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I. INTRODUCTION 

The presence of large scale magnetic fields in galaxies, galaxy clusters and in intergalactic 
space [1] is a mystery and several theoretical models have been presented to explain the 
origin of these fields [2-5]. Most of these works deal with the dynamo theory of single 
fluid magnetohydrodynamics (MHD). But the set of MHD equations assumes that some 
magnetic field is already present in the system. Therefore, these models actually investigate 
the amplification of the existing weak magnetic field and can not explain the generation of 
the 'seed' field in true sense. 

Long ago [6] , Biermann presented a mechanism for the generation of stellar magnetic fields 
which is not based on MHD. He proposed that the electrons faster motion compared to ions 
can produce electric field in a rotating star which is not curl free due to non-parallel density 
and temperature gradients and hence the magnetic field is produced. The ions were assumed 
to be stationary in this work. The Biermann battery and electron diffusion processes were 
also investigated to explain the generation of 'seed' magnetic fields in galaxies [7] . 
It is very interesting that the large magnetic fields of the order of kilo and mega Gauss were 
observed in classical laser- induced plasmas many decades ago [8-9]. These observations 
indicate that the dynamics of initially unmagnctized nonuniform plasmas can generate 
magnetic fields. The idea of magnetic field generation by plasma dynamics is very attractive 
and a huge amount of research work in this direction has already been appeared in literature. 
Based on Biermann battery effect, a single fluid plasma model called electron magneto- 
hydrodynamics (EMHD) was presented to explain the magnetic field generation in laser 
plasmas [10, 11]. But it does not require the rotation of the system to produce magnetic 
fields. In EMHD, ions are assumed to be stationary and electrons arc treated to be 
inertialess. Both the fluctuating [12, 13] and steadily growing magnetic fields [8, 9] have 
been theoretically produced using EMHD models. The so called magnetic electron drift 
(MEDV) mode was discovered [12] using EMHD equations. The MEDV mode is believed to 
be a pure transverse low frequency wave of an unmagnetized inhomogeneous plasma. But a 
critical analysis of MEDV mode shows that it should contain a contribution from electron 
density perturbation as well. 

Therefore, a new mode which is partially transverse and partially longitudinal has been 
proposed to be a normal mode of pure electron plasmas which can exist only in a very 
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narrow range of parameters. This mode can couple with the ion acoustic wave which also 
becomes electromagnetic under certain conditions in a non- uniform plasma [14]. Similarly 
the model equation widely used for the generation of steadily growing magnetic field [12, 13] 
is not fiaw-less. The field grows on ion time scale while ions are assumed to be stationary. 
Recently [15], the same model equation containing electron baroclinic term has been used 
to estimate the magnetic field produced in a laser plasma. 

Some weaknesses and contradictions in the approximations and assumptions used in EMHD 
models for magnetic field gereration have already been pointed out [16]. The advantage 
of EMHD model is that it is very simple. Since it is still being used by many authors, 
therefore it seems important to discuss at least the two cases of fluctuating and steadily 
growing magnetic fields which are believed to be generated by EMHD. The EMHD models 
for the generation of magnetic fields are critically discussed in the next section. 
The magnitudes of magnetic fields on galactic scale [7] as well as on laser-plasma scale [8, 9] 
were estimated by assuming the non-parallel electron density and temperature gradients to 
be constant using EMHD equations. In these models the gradients of electron temperature 
and density were assumed to be one-dimensional and the produced magnetic fields had only 
one component. 

A few years ago [17], a theoretical model was presented to show the generation of three- 
dimensional magnetic field by baroclinic vectors of electrons and ions. However, in this 
investigation some constant magnetic field was assumed to be present already. A stationary 
solution of these equations was presented by Mahajan and Yoshida [18] in the form of 
double Beltrami field. 

Later, the model presented in Ref. [17] was modified to explain the creation of all the three 
components of 'seed' magnetic field vector B from t=0 due to externally given forms of 
barochnic vectors [19]. Here one does not need to assume some static magnetic field to 
be present in the system. However, the form of the solution was sinusoidal along one axis 
which is not physical in general. 

The EMHD may have some applications in other areas but for the magnetic field generation 
on slow time scale, the ion dynamics can not be neglected. Therefore, it is necessary to 
study the 'seed' magnetic field generation by using two fluid model. On the other hand 
the generated magnetic field vector should not have necessarily only one component for the 
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sake of generality. 

Our aim is to find out an analytical solution of the set of two fluid equations such that the 
cross products of plasma density and temperature gradients of electrons and ions become 
the source terms in the electron and ion equations of motion. We assume that the plasma 
has been produced in a non-equilibrium state and it evolves with time generating the "seed" 
magnetic field and fiow. 

The present investigation is very different from the previous work [19] because the density 
gradient scale length is assumed to have an exponential form in xy-plane. In most of the 
analytical studies, the density is assumed to have exponential form along one axis. But we 

want to obtain a two-dimensional solution, therefore the density is assumed to be a function 
of (x,y) coordinates. 

We have chosen special profiles of gradients of density (Vn) electron temperature (VTe) 
and ion temperature (VTj) to obtain an analytical solution. Different profiles of density and 
temperatures can be considered but then the numerical simulation will be needed. In our 
formulism all the nonlinear terms vanish and ultimately we obtain the two linear equations 
where the terms VV' x VTj (i = e, i) with -0 = Inn (where n is normalized density) become 
the source terms for magnetic field. 

The details of the model are discussed in section III. Since the present model contains a very 
complex system of highly nonlinear equations, therefore to find out an analytical solution 
one has to use some assumptions and approximations. The main focus is to justify the 
physical idea that the two fiuid plasma with density gradient like an exponential function 
in xy-plane and constant gradients of electron and ion temperatures along y and z-axes can 
generate the 'seed' magnetic fields and flow. This model can explain the magnetic flelds 
produced in laser-induced plasmas. In our opinion, the numerical simulation of two-fluid 
equations is very important to study the 'seed' held generation. In simulation one can use 
many different proflles of density and temperatures. 

II. CONTRADICTORY RESULTS OF EMHD 

Here we briefly point out the contradictory results of EMHD models used for magnetic fleld 
generation. First we discuss a theoretical model based on EMHD for the generation of 
fluctuating magnetic flelds proposed several years ago [12]. The mode-discovered through 
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EMHD theory was named as the magnetic electron drift vortex (MEDV) mode. A great deal 
of research work on this mode has been carried out. 

The critical analysis of MEDV mode and some new theoretical results have been published 
recently on the fluctuating magnetic fields [14]. There seems to be a need to briefly clarify 
here the physical situation to lay down the basis of our theoretical model presented in the 
next section. In the theory of MEDV mode, the ions are assumed to be stationary but the 
electron inertial effects are included. The hnear description of the MEDV mode is presented 
very briefly as follows. 
Electron equation of motion is. 



meUodtVei = -enoEi - Vpei 



(1) 



where subscripts one (1) and naught (0) denote the linearly perturbed and equilibrium 
quantities, respectively. In the limit oUpi « ou « oUpe (where uj is the frequency of the wave 
and Upj — { ^^^.^ ^ is the plasma oscillation frequency of the jth species while j= e here), 
the displacement current is ignored and Maxwell's equation yields. 



47r 47r 
V X Bi = — (Ji) = — (-enoVei) 
c c 



(2) 



Since V.Ji = 0, therefore according to (2), the density perturbation is neglected and we find 
Pel = noTf-i. For Tgi, the electron energy equation becomes. 



3 3 

-nodfT^i + -no(vei.V)Teo = -poV.Vgi 



(3) 



Assuming, Vno = x I ^ I , 



1 dno 
no dx 



, El = £^ix, k = kyy and Bi = Biz the linear 



dispersion relation for MEDV mode turns out to be 



2 2 , x2 2 I 2 



(4) 



where Cn = , ,%l2 , Ag = — is the electron coUision-less skin depth and vxe = — is 
the electron thermal speed. If temperature gradient is assumed to be anti-parallel to the 
density gradient in laser plasma with VTq = x |^| and kt = ■, then (4) is modified 

(5) 



as, 



Co 



1? p 
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and these magnetic perturbations become unstable if the condition 

2 

-K,n< f^T (6) 

holds. Note that the local approximation requires kt « ky. It has been assumed that 
V.Ei = and V.Vei 7^ in the description of MEDV mode along with Upi « u. 
Equation (1) indicates that the term Vpei = V(noTei) = T^iVtiq + noVTgi will produce a 
linear term with V replaced by k and hence Ei can have a longitudinal component with 
V.Ei 7^ 0. Thus the mode can not be a pure transverse mode. 

Moreover, the linear theory has been applied under the local approximation therefore the 
term, {jr^ "^rj^y ^^"^ be closer to c^/c^ where Cg — {Te/rriiY/'^ is the ion acoustic speed while 
Co < 1 always holds. Using the laser plasma parameters [9] Te — lOOey and no ~ 10^° cm"^, 
one obtains ui < uipi contrary to initial assumption of stationary ions for uipi « ui. 
Recently [14], it has been shown that if compressibility effects arc also taken into account 
then one obtains a new partially transverse and partially longitudinal normal mode of a 
nonuniform pure electron plasma as, 



where H,^[{l + fA^^,'} a - ^l/k^] and a = (1 + Xlk^). 

It is important to note that ions can be assumed to be stationary in the limit ^ — > 0. 
Since for hydrogen plasma ^ ~ 10~^, therefore (7) is valid for me/rrii < }?£,g^ky, i^'i,/ky and 
« ul^. 

But the important point is that the term ^'l^g'^n can be closer to c^/c^ and lesser than u)"^ 
therefore the dynamics of ions should not be ignored. 

It is also interesting to mention here that ion acoustic wave (lAW) has always been treated 
as a low frequency electrostatic mode. The reason is that in the limit ^ — > 0, the inertia- less 
electrons are assumed to follow the Boltzmann density distribution in the electrostatic field 
E = as. 



~ e-^^l^^ (8) 



In an inhomogeneous plasma we may have ^ < i\n/ky and in this case electron inertia 
should not be neglected. Then for ^ < Xj^^ky, the lAW follows the dispersion relation [24], 



^2 ^ c2i;2^I -'n, -yj /gx 
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where b — (l + Xj^^ky) . Hence inhomogeneous plasmas can have a low frequency electro- 
magnetic wave on ion time scale. 

When the electron temperature perturbation effect is taken into account, then modes de- 
scribed in (7) and (9) will couple to produce a partially longitudinal and partially transverse 
wave with the dispersion relation [24], 



If Vpeo = is used as the steady state condition the above equation yields a basic low 



frequency electromagnetic wave of inhomogeneous unmagnetized plasmas with the dispersion 
relation, 

R r . . . . / _i 



3//o 

This wave has not been studied in plasmas so far. In our opinion it can play very important 
role in the generation of magnetic fluctuations in unmagnetized plasmas due to several 
linear and nonlinear mechanisms. In a pure electron plasma where ions are assumed to be 

stationary, (10) reduces to 

-^'=^yekl)« (12) 

But in our point of view the electron plasma wave frequency in (12) is near ion acoustic 
wave frequency c^gky and hence it couples with it. 

Now we look at EMHD theory for the generation of 'seed' magnetic field which is steadily 
growing. Again ions are assumed to be stationary in the time scale r << u}~l. In addition 
to this the electron inertia is also neglected assuming « r . Then electron equation of 
motion becomes, 

~ -eE - ^ (13) 
n 

The faraday law is, 

dtB = -cV X E (14) 
If it is assumed that Vrio = x and VTq = y |^^|, then (13) and (14) yield. 
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where L„ — k^^ and Lt — Kj} are constants. 

Equation (15) is integrated from r = to r = ^ to have [9 - 11], 




e \LnLT 




(16) 



This is a well-known equation in laser-plasma literature. Assuming Tg ~ lOOey, 
no ~ 10^° cm~^, Ln ^ Lt ^ 0.005cm, one obtains Cs ~ 3 x lO'^cm/Sec and hence 
\B\ - 0.6 X 10^ Gauss [9]. Note that r = ^ ~ 1.66 x 10-^° Sec and oopi ~ 1.3 x lO^had/Sec 
while the laser pulse duration is of the order of a nano second. Thus we have uj~^ « r 
contrary to initial assumption of stationary ions for r << u)^- An equation similar to (16) 
has also been used to estimate the 'seed' magnetic field generated by ionized clump of a 
galactic cloud [7] . The density gradient of the cloud has been assumed to have exponential 
form. 

The brief overview of EMHD models shows clearly that the theoretical models for both the 
fluctuating and steadily growing magnetic flelds suffer from serious contradictions. Since 
the EMHD is still being used [15] for estimating magnetic fields produced in laser plasmas, 
therefore the weaknesses and contradictions have been elaborated here again. 
Since the plasmas have generally exponential density profiles, therefore there is a need to 
find out an exact 2-D solution of the set of two fiuid equations assuming exponential type 
density structure in a plane. It may also be mentioned here that in many tokamak plasmas, 
the density falls almost exponentially near the walls which gives rise to drift waves. 

III. EXACT SOLUTION OF 2-FLUID EQS. 

Our aim is to search for an exact analytical solution of the set of highly nonlinear partial 
differential equations of electron-ion plasma to show how the system from a non-equilibrium 
state can evolve in time generating the 'seed' magnetic field. 

In our opinion the same physical mechanism is applicable at both astrophysical and 
laboratory scales. Biermann [6] gave the pioneering idea that the electron baroclinic vector 
(Vrie X VTg) can generate magnetic fields in rotating stars. We just modify it a little by 
proposing that the 'seed' magnetic field is a macroscopic phenomenon and it is generated 
on longer spatial and temporal scales. Hence the ion dynamics can not be neglected. 
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Therefore, the ion baro-chnic vector (Vn^ x VTj) is also crucial to be considered. However, 
the quasi-neutrality approximation in slow time scale is a valid approximation, therefore we 
use rij ~ TT-e = n. In thermal equilibrium, the source of magnetic field generation disappears. 
For an analytical solution of a complex set of equations, we have to use some assumptions 
and approximations. The exact solution presented here contains exponential type density 
profile in xy-plane which is the main deviation from the previous models [17, 19]. It is 
important to note that the nonlinear terms do not vanish if we assume exponential density 
fall or rise along both the axes x and y as has been discussed in section II. 
Therefore, we have to choose a very special exponential function in xy-plane for density 
which reduces the nonlinear equations into two linear equations. The detailed mathematical 
model is presented here. 



The electrons are assumed to be inertialess in the limit \dt\ « cOpe, c|V| where cOpe — 
( ^^m°^^ ) ^ electron plasma frequency and c is the speed of light. We define four scalar 

fields ip, u, X and h such that the ion velocity Vj and magnetic field are defined, respectively, 
as [17, 19], 

Vi = {V(p X z + uz)f{t) = {dy(p, -d^(p, u)f (17) 
B = (Vx X z + h7)f{t) = {dyx. -d.X, h)f (18) 

All these scalar fields are functions of (x, y) coordinates and f is a function of time. 
We further assume dtUj = and V.Vj = which requires 

VV'.v, = {(^,V} = (19) 

where ijj = Inn, ne — ni = n and {cp, ijj} = dy^pd^ip — dx^pdyijj. Here n = "^^'^^ and A^o is an 

arbitrary number used to normalize n. 

The displacement current is neglected and hence we obtain, 

V, - (20) 



47re n 

Let E = — V$ — ^dtA where $ is electrostatic potential different from (p. Since B = at 
t=0, therefore we do not normalize the equations. This point has been explained in detail 
in Ref. [19]. The curls of momentum equations of electrons and ions yield, respectively, 

c V X B 



B X 



4:716 n 



^(W X VTe) (21) 
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and 

dt{aB + V X Vi) - V X [a{vi x B) + x (V x v^)] 

= —W X Ti) (22) 

where a — If the conditions 

{^,x}^W,u}^{h,^}^0 (23) 

satisfy along with 

{ VV, = (24) 

then all the nonlinear terms of (21) and (22) vanish and they reduce, respectively, to simpler 
equations 

dtB = --(VV^ x VTe) (25) 
e 

and 

dt(aB + V X Vi) = — (VtA x VTA (26) 
rrii 

where Te ^ % and right hand sides of (25) and (26) are the source terms for generating 
magnetic field and plasma flow. Let us assume that B is related with plasma vorticity 
through the following equation, 

B = a (V X Vi) (27) 
where a is a constant. Then (26) becomes 

(a + a-^) dtB = — (V^^ x VTA (28) 
rrii 

Now we discuss an important point of the present theoretical model. In the previous works 
it was assumed that the field (p satisfies the Poisson equation, 

VV = -A<^ (29) 

where A is a constant and < A holds. The forms of V' and Tj were chosen as = 

il)Qe^^^cosiJL2y and Tj = {Tooj +T'Qj{y — z)} f {t) where ■?/'o, /^i, /^2, Tooj and Tq^- were constants. 
Wc, here, want to find out a 2-D solution in the exponential form without assuming the 
density gradient to be constant. For this purpose the assumption (29) is modified as 

VV = Ay? (30) 
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and we assume < A. The form of can be chosen hke, 

^{x,y) = Aie^^^+'^y') + A2e^^''-''y^ = + = Inn (31) 

where n = ""^^"^ is dimensionless and A^o is some constant density. Here Ai, jji, A2 are 
constants and 

We may choose < for simphcity. Let the temperatures be only functions of space in 
yz-plane as, 

Toj {y,z) = Tooj + Tqj {y - z) (33) 

Then the barochnic vectors in (25) and (26) become constant with respect to time. These 
equations can be integrated from t=0 to r and one obtains, 

B = -^(W X V7;)t (34) 

and 

B = ^ _ ( W X VTOr (35) 
mi{a-\-a ^) 

These equations relate T'^q and T^'q as. 

Equations (31) and (33) yield 

{V^ X VT,) = -To,, {dy^, -d.ij, -d^^P) (37) 
where dyi}) = v {ipi — ■^2) and d^ip = A^V'- 

We arc free to use any of the equations, (34) or (35) to estimate B. Let us choose (34) and 
use (18) for f=l to find out following relations, 

X = (XoV')t (38) 

and 

h = -(/io'0)t (39) 
where Xo = {f~e^^ ^0 ~ '"^o- Equation (17) for f=l along with (30) gives, 

V X Vj = {dyu, -d^u, -Xip) (40) 
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Then (27) yields, 

u = (lioV')T (41) 

and 

ip = M)t (42) 

where uq = ^ and v^o = (ix) ^o- 

Three dimensional magnetic field B and Vj can be expressed explicitly for Ai ^ A2 as, 



B 



and 



V,; 



-00 



(43) 



i/(^o(Ae('^^+'^^) - Ase^'^^-'^^')) 
-/i(/?o(>lie(''''+''«') + ^126^'^'^-''^)) 
-uo{Ale^^''^+''y^ + A2e^^''-''y^) 



(44) 



Hence all the scalar fields x, h, u and </? become functions of (x,y) through which is 
externally given. The complicated nonlinear terms of two-fiuid equations vanish and we 
obtain two simple and beautiful linear equations (34) and (35). These equations show that 
the terms of electrons and ions {Vip x VT^) (for j = e, i) become the source for the 'seed' 
magnetic field and plasma fiow Vj. 

IV. GENERAL APPLICATIONS 

It has been shown that the forms of ip{x, y) and Toj{y, z) given in equations (31) and (33), re- 
spectively, reduce the set of nonlinear two fluid partial differential equations into two simpler 
linear equations (34) and (35) under certain conditions mentioned in the previous section. 
This theoretical model shows that the non-parallel density and temperature gradients can 
create magnetic fields and fiows in initially unmagnetized plasmas. 
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Our aim is to apply this model to a system which has exponential-type of density profile as 
for example in the case of a laser plasma. But the chosen form of ip{x, y) in (31) with the 
definition ip — Inn gives, 

where Aq — Ai — A2 has been assumed. 

It looks as if density n(^j:^y) has a profile like double exponential in (x,y) plane. Such a steep 
density variation is not interesting for physical applications, in general. We show here that 
the density variation can become very similar to exponential form in (x,y) plane by choosing 
suitable values of the constants A^o, A* ^ along with Aq. Let us consider an inhomogeneous 
plasma rectangle in (x,y) plane with four corner points (0,0), (x^, 0), (0, y^) and {xm,ym) 
where Xm and ym are the maximum lengths of the system along x and y axis, respectively. 
Then choose the constants in such a way that the density n(^x,y) at {xm,ym) will be almost 
e-times (or a little larger) than the value at (0,0), while the density at (a;^,0) and (0,1/^) 
will be somewhat lesser than e-times the density at (0,0). 

Such a density function is acceptable physically. For example, in previous EMHD model, the 

X 

density was chosen to be an exponential function only along x-axis as n(^x) — UqC^ where 
Ln is the density scale length and no is the magnitude of density at x=0 [9]. In our case, n 
depends upon two coordinates x and y. It's profile depends upon the values of the constants. 
Therefore this theoretical model can be applicable to many inhomogeneous plasma systems. 
Note that ■?/' = Inn and if < < 1, < + z/y < 1 and < fix — uy < 1, then ip 
can have values which give e'^ of the order of e^^\ but not exactly e(l) ~ 2.7 at all points 
because ip changes with x and y. 

In the next section we shall apply the model to laser plasmas as an example and our point 
of view will become clearer. The values of the constants will be chosen to show how it works 

for relatively smooth density profiles in (x,y)-plane. 

It seems important to point out that any one of %bi and 1^2 in (31) should not be much 
smaller than the other while choosing the constants. If one of them is negligibly small then 
the solution becomes one-dimensional. 

V. LASER PLASMA 
Here we apply our theoretical model to estimate magnetic field B and the plasma fiow Vj 
in a non-uniform classical laser plasma. Consider a finite plasma rectangle with four corner 
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points (0, 0), {xm, 0), (0, y^) and {x^, Um) in (x, y) plane as has been mentioned in previous 
section. We may assume /ix and uy to vary in this finite plasma as, 

//x : — > 0.5 = fj,Xm (46a) 

and 

i/y : ^ 0.7 = (466) 

Then at (0, 0), we have, n(o,o) = ^^(^ is chosen such that n 1 or n 1 because 
density n(^a:,y) should be neither zero nor negative. Therefore, we choose = 3 which 

gives, ^(0,0) — 1-1 and due to (45) we find Aq ~ 0.55. If density is of the order of 10^° cm~^, 
then we may assume A^o — 10^° cm~^ and hence 7^(0,0) = 3 x 10^° cm~^. Or we may assume 
A^o — lO^^cm"^ and hence 7^(0,0) = 3 x lO^^cm"^ while n(xm,ym.) ^^i^^ ^^rn out to be nearly 
10^°cm~^. In laser-plasmas, = 50 x 10~^ m — Lt was assumed in estimating \B\ [8, 
9] (where L„ and Lt are scale lengths of density and electron temperature along x and y 
co-ordinates, respectively). In these studies, the Vn was along x-axis and VTg was along 
y-axis only. Then L„ ~ Lt was also assumed. For the sake of generality we do not assume 
II — u. Instead let 1/ — 1.5/1. In this case we estimate n, and B at 4-corner points of the 
plasma rectangle as follows: 

n(o,o) = 3, ^'(0,0) ^ 1-1 
B(o,o) = (0, -1.1, -1.1)9.9 X IQ^Gauss (47a) 
n(x„,o) = 6 : ^'(^„,o) ^ 1-79 
B(^^,o) = (0, -1.79, -1.79)9.9 x lO^Gauss (476) 
nio,y^) = 3.98; "^(y^fl) = 1.38 
^{o,ym) = (1-24, -1.38, -1.38)9.9 x lO^Gauss (47c) 

^(a..,.^.)- 9-48; *(.^,^^)c^ 2.25 
^{xm,ym) - (2.73, -2.25, -2.25)9.9 x lO^Gauss (47d) 
These values of |B| are almost in agreement with the observation [9]. Then we can express, 

Vi(^,j,) = 6 X 10^(-0.46, 0.3, -l)^(^,j,)cm/sec (48) 
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If we look at the values of density, we notice that density n(^x,y) {^m, Vm) is 

H^m,ym) - {^(0,0)} (3. 16) 

which is little larger than e-times the density at (0, 0). Therefore, the density has not 
been chosen as a double exponential function. We are mainly interested in the order of 
magnitudes of |B| and |vj| which mainly depends upon our choice of constants. 

V. DISCUSSION 

It is important to note that the 'seed' magnetic field generation can not be explained on 
the basis of magnetohydrodynamics (MHD). There must be a source like thermal energy 
which converts into magnetic energy. Biermann [6] proposed that the electron baroclinic 
vector (Vrie x VTe) can generate the magnetic fields in rotating stars. Then based on this 
idea, a very simple model; the electron magnetohydrodynamics (EMHD) was presented to 
explain the generation of magnetic field in laser-induced plasmas. Later on, the magnetic 
electron drift vortex (MEDV) mode was discovered using EMHD. This was beheved to be a 
low frequency pure transverse normal mode of electron plasma [12, 13]. The EMHD is also 
not a convincing theoretical model for the generation of 'seed' magnetic field. The MEDV 
mode description also suffers from contradictions. 

On the other hand, the two-fiuid model is too complicated. The numerical simulation of these 
equations is a complex problem. But the analytical 3-D solution is also not straightforward. 
However, a two-dimensional solution has been presented using physical and consistent as- 
sumptions and approximations. It is important to note that the present model is different 
from the previous works because it shown that 

1. ion dynamics play a crucial role 

2. the baroclinic vectors generate not only the magnetic field but plasma fiow as well 

An exact 2-D solution of the two-fiuid equations was also found a few years ago [19], but 
it had a serious weakness. The density gradient was assumed to follow sinusoidal behavior 
contrary to common observations. Since, it was the first effort to get an exact solution of 
the two fiuid equations, therefore it was presented for the interest of researchers working in 
the field. 
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The present 2-D exact solution of two-fluid equations is in the form of exponential function 
of X and y coordinates. This structure of density gradient is more physical. Since all flelds 
become linear function of ijj, therefore all fields have the similar spatial structure. This 
solution is applicable to both astrophysical and laser-induced plasmas in our opinion. 
The present investigation suggests that instead of EMHD, the numerical simulation of two- 
fluid equations will be very useful for understanding the mechanism for the generation of 
'seed' magnetic fleld in different systems with different proflles of density and temperatures. 
For analytical solution, we have to choose very special forms of density and temperature 
gradients. This theoretical model can be very useful for further studies in astrophysical and 
laser plasmas. 
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